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We calculate supercurrent across a two-dimensional topological insulator subjected to an external magnetic 
field. When the edge states of a narrow two-dimensional topological insulator are hybridized, an external mag¬ 
netic field can close the hybridization gap, thus driving a quantum phase transition from insulator to semimetal 
states of the topological insulator. We find a sign reversal of the supercurrent at the quantum phase transition 
revealing intrinsic properties of topological insulators via Josephson effect. 
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I. INTRODUCTION 

Two-dimensional topological insulators (TIs) in the pres¬ 
ence of time reversal symmetry are insulators in bulk while 
they possess gapless metallic edge states [1, 2]. These edge 
states are helical, i.e. they support counter-propagating modes 
with opposite spin projections, and are protected against elas¬ 
tic backscattering. One interesting phenomenon that may 
arise due to the rigid spin-momentum locking is topological 
superconductivity [3-6]. The topological superconductivity 
can play an important role in the realization of non-Abelian 
statistics needed for topological quantum computation [2, 7]. 
To establish superconducting edge states in TIs, one can prox- 
imitize these states with a superconductor which results in the 
forming of Cooper pair wave functions in the TI (reviews on 
the proximity effect can be found in Refs. 8 and 9, for in¬ 
stance). 

The edge states of a two-dimensional TI were exper¬ 
imentally investigated in two different material classes 
(HgTe/HgCdTe and InAs/GaSb) through imaging magnetic 
fields produced by corresponding edge currents and determin¬ 
ing an effective edge resistance using a scanning supercon¬ 
ducting quantum interference device (SQUID) [10, 11]. In or¬ 
der to reveal the topological superconductivity in these edge 
channels, the current-voltage characteristics and Fraunhofer 
interference patterns in superconductor - two-dimensional TI 
- superconductor junctions were recently studied both theo¬ 
retically [12-17] and experimentally [18-21]. In particular, 
the experimental study of the Fraunhofer interference patterns 
may allow one to argue for an edge-dominated contribution to 
the Josephson current although not conclusively [20, 21]. 

An important practical challenge in Refs. 20 and 21 might 
be the penetration depth of helical edge states into the bulk 
material. According to these experiments, the penetration 
depth could be comparable to the width of the TI strip. There¬ 
fore, the edge states on the opposite sides of the strip may hy¬ 
bridize and produce a gap in the spectrum of edge states [22]. 
However, the externally applied magnetic field normal to the 
TI plane can close this gap and result in a phase transition 
from insulator to semimetal [22-25]. Here, we demonstrate 
that the latter apparent disadvantage can turn to give evidence 
for the existence of edge-mode superconductivity in a point 
contact Josephson configuration as sketched in Fig. 1. We 
show that, by tuning the chemical potential, such a hybridiza¬ 


tion between the edges in the presence of external magnetic 
field results in a phase transition and reversal of a supercur¬ 
rent flowing across the edge channels. 

In what follows, we first describe the model of the Joseph¬ 
son junction, discuss the properties of the two dimensional 
TI band structure, and then illustrate how the interplay of the 
magnetic field and hybridization between the edge states af¬ 
fects the supercurrent. 


II. MODEL OF THE JOSEPHSON JUNCTION 

The Josephson point contact configuration we consider in 
this paper is schematically shown in Fig. 1. We choose the 
plane of the two-dimensional (2D) TI to reside in the z = 0 
plane, the top and bottom edges of the TI to be parallel with 



FIG. 1. (Color online). Schematics of the point contact superconduc¬ 
tor (SC) - 2D topological insulator - superconductor junction. The 
topological insulator strip has a width of d and the junction plane 
resides in the z = 0 plane so that the superconducting contacts are 
located at r-p/B = ±d/2, 0). The external magnetic field 

H = (0, 0, H) is oriented along z, normal to the plane of the TI. 
Top and bottom edge states (shown by parallel arrow lines) are hy¬ 
bridized while left and right edge channels are strongly gapped com¬ 
pared to the region of contact. The top and bottom superconducting 
electrodes have different macroscopic phases marked by fx and (j>B, 
respectively. The typical tunneling process of Cooper pairs between 
the superconductors through the edge channels is also illustrated. 
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the X axis. The two superconducting electrodes are connected 
to the TI via top and bottom edges at = (xt, +dl2, 0) and 
Vb = {xb, —d/2, 0), respectively. An external magnetic field 
is applied normal to the junction plane along the z direction, 
H = (0, 0, H) (we will assume H > 0 henceforth). We adopt 
the Landau gauge for the vector potential A = {—yH, 0, 0) 
throughout our calculations. In what follows, we will be using 
h = 1 units. Therefore, the system shown in Fig. 1 can be 
described by the following Hamiltonian: 

H = Hsc + HtI + HTun, (1) 


where the first term is the usual BCS Hamiltonian describing 
superconducting leads 


dfsc = 


j 




— (-fVr + eKjcf 
Zm 





Aj(r)T>]_^(r)$]_^(r) +h.c. 



where ^ (r) is the electron creation operator in an SC lead, 
j = T/B labels top/bottom SC lead, a =t, i stands for the 
spin quantum number, p, = pjp/2m and m are chemical po¬ 
tential and electron mass, and (r) is the superconducting 
order parameter. The integration is performed over volume 
Vt/b of the top/bottom SC. We consider magnetic field to be 
smaller than the first critical field of the SC so that orbital ef¬ 
fects are weak and also neglect the Zeeman effect in the SC. 

The second term in Eq. (1) stands for the Hamiltonian of 
the TI strip. 


Hti = 




dx'^jc, 


(x) 


+ tTjj,6aa' - 


33'^ 
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( 2 ;), 


( 3 ) 


where is the electron creation operator in the 

top/bottom edge states, v is the Fermi velocity, characteriz¬ 
ing the edge dispersion, t is the tunneling amplitude between 
top and bottom edge states ( which ensures electron spin and 
momentum conservations and we choose f > 0 for concrete¬ 
ness), CT* and T®, i = x,y, z, are the Pauli matrices acting on 
the spin and the top/bottom edge of the TI pseudospin degrees 
of freedom, respectively. The position of chemical potential 
in the TI is measured from the charge neutrality point and is 
defined by Sy,. We have also defined a magnetic wave-vector 

kh = d/2i^ - gyBH/2v, (4) 


containing contributions from the Aharonov-Bohm phase gra¬ 
dient (first term) and Zeeman coupling (second term), in 
which i = \JcjeH is the magnetic length, g and p b are the 
electron g-factor and Bohr magneton, respectively. 

The third term in Eq. (1) describes the tunneling of elec¬ 
trons between superconducting leads and edge states of TI, 


^Tun — ^ ^ / dx 


-d/2, 0)T'B,a(a;) 


-h to{x)'^\,^^{x)^T,a{x, d/2, 0) -P h.c. 


( 5 ) 




FIG. 2. (Color online). Band structure, Eq. 6, of the 2D TI for dif¬ 
ferent values of the magnetic energy eh and the tunneling amplitude 
t. Arrows indicate the spin projection of an electron in the corre¬ 
sponding subband, (a) The hybridization between the helical edge 
states with the same spin opens a gap of 2t at zero momentum. The 
subbands are doubly spin degenerate in the absence of an external 
magnetic field, (b) The external magnetic field removes the spin de¬ 
generacy and closes the gap. (c) Critical point eh = t where two 
bands touch at fc = 0. (d) Subbands with opposite spins cross at two 
points when eh > t. 


where tQ{x) is the tunneling amplitude between the supercon¬ 
ductor and the edge states. In our calculations we consider the 
geometry depicted in Fig. 1 where the electron spin is con¬ 
served in the tunneling process between superconductor and 
TI whereas its momentum is not. 


III. SPECTRUM OF HYBRIDIZED EDGES STATES IN 
MAGNETIC FIELD 


Now we discuss the properties of the TI edge band structure 
in the absence of superconducting leads. We shall return to 
the calculation of Josephson current later. Diagonalizing the 
Hamiltonian (3) we obtain the edge states band dispersion, 

Es,s' (k) = SEh + s'\/ -I- ^2^ ( 5 ) 

where s, s' = ±1, k denotes the momentum along the edge, 
and Eh = v\kh \ is the “magnetic energy”. Figure 2 illustrates 
the characteristics of the band structure for various values of 
Eu/t at the neutrality point 5y = 0. In the absence of an 
external magnetic field {eh — 0) the spectrum consists of sin¬ 
gle valence and conduction bands split by a gap of size 2t at 
fc = 0, as shown in Fig. 2a. In this regime, the system is an 
insulator and each band is doubly degenerate in spin space. 
However, a small magnetic field Eh < t lifts the spin degen¬ 
eracy and splits both conduction and valence bands by 2eh 
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at fc = 0 which is seen in Fig. 2b. At a critical field value, 
namely eh = t, the lowest conduction and highest valence 
bands touch at a single point, k = 0. When Eh > t the sys¬ 
tem is gapless and the low energy part of the spectrum has two 
crossing points at ±fco, where ko — \//v, see Fig. 
2c,d. Thus, the critical value of the magnetic energy, eh = t, 
defines a quantum phase transition between the semimetallic 
and insulating states of the TI narrow strip [22, 23]. Below, we 
show that the insulator to semimetal phase transition can have 
nontrivial signatures in the Josephson current flowing through 
such a system. 


IV. EVALUATION OF THE JOSEPHSON CURRENT 

We assume that the tunneling amplitude to between the TI 
edge states and superconducting electrodes is smaller than the 
one between the edge channels. This assumption allows us 
to neglect the proximity induced superconducting gap in the 
TI. The Josephson current J in the lowest order in to can be 
expressed as 


J = 2er^ ^I|(fa:Jm|FB_aa'(w„; 


2 = 1 


Xl,X2)to{x2) 


X Gtbm' (-w„; Xo, X 2)to{xo)Fl.^^,^{uin\Xo, Xi) 


X to{x4)GTB.a{,OJn]X4,Xi)toixi) 


(7) 


where = (2n + l)7rT is the fermionic Matsubara fre¬ 
quency, n G Z, and T is the temperature. More details on the 
derivation can be found in Appendix. To derive Eq. (7) we 
have assumed that the coupling of superconducting leads via 
edge states along the perimeter of the TI strip is strongly sup¬ 
pressed in comparison to the direct tunneling across the TI. 
Such a condition can be experimentally achieved, for exam¬ 
ple, by inducing an insulating gap in the left/right edge states 
stronger than that of the region in the vicinity of the contacts. 

The anomalous Green function inside the superconduc¬ 
tor Fj_aa'(T, r', r, r') = (here the 

operator „ (t, r) is in the Heisenberg representation) in 
the frequency representation and quasiclassical approximation 
reads. 


77,>c<'(wn;r,r') = 


dp 

(27r)3 + 


where fp = p^/2to — /i, the absolute value of the order pa¬ 
rameter |A| is assumed to be fixed in both superconductors, 
and (j)T/B is the macroscopic phase of top/bottom lead. 

The hybridization of the helical edge states is de¬ 
scribed by the electron Green function GTB,a{x, r'; x, x') = 
— {TT^T,a{j,x)^\^{T',x')), which can be obtained from 
Eq. (3) in the frequency representation as follows. 


GTB.t/f 2 ;, x') 


dk/{2n) 

(iujn -f d/r ± eh)^ — — F ' 



FIG. 3. (Color online). Normalized critical supercurrent 
JcisH, Sfi)/Jc{0,0) through the point contact configuration pro¬ 
posed where |A|/t = 6. (a). Critical current as a function of 
magnetic energy normalized hy the tunneling amplitude at differ¬ 
ent values of chemical potential: S^i/t £ [0,0.6,1,1.4]. (h) Crit¬ 
ical current as a function of chemical potential dp normalized by 
the tunneling amplitude t for different values of magnetic energy: 
SH/te [0,0.5,1,1.5]. 


It is worth noting that Eq. (7) governs the direct tunnel¬ 
ing process of two electrons with opposite spins across the 
hybridized edge states. We adopt a tunneling barrier model 
where the hopping amplitude to{x) = to6{x — Xt/b) is 
considered for top/bottom superconducting contacts in which 
xt,b G [—L/2,L/2]. This assumption allows us to inte¬ 
grate the Josephson current, Eq. (10), over the lengths of 
top and bottom edge states involved in the supercurrent trans- 
port: (J) = J_^^^dxTdxBJ/■ This double integral can 
be rewritten in terms of a double integral over the relative and 
center of mass coordinates. 

Substituting the Green functions from Eqs. (8) and (9) into 
Eq. (7) we obtain, 

J = 4e [7r^jzTii/sc|A|tto]%in((/)B - cjir) 

_^ 1 g—|2 :b— 

^ n ^; 2_|_|^|2 _ jg ^)2 _|_ ^ 2 Jl /2 

( 10 ) 

where bsc = mpp/^ir^ is the density of states at the Eermi 
level per spin in the normal state of the superconductor leads, 
and bti = 1 /2Trv is the density of states per spin and TI edge. 
We are interested in the zero temperature limit, T = 0, and 
transform the sum over frequencies into an integral, 
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/ ir- emphasize that the integral in (J) behaves qualita¬ 
tively different in the presence and absence of band crossing 
points in the spectrum, see Fig. 2. Indeed, we obtain a simple 
analytical expression for the Josephson current in the narrow 
junction limit L v/t of the form 

{J) = Jc{eH,Sfi)sin{(j)B - (jir), (H) 

where the critical cutTent across the contact at the charge neu¬ 
trality point, Sjj. = 0, is given by, 

Jc{£h, 0) = e| A|f [TT^J^TiJ^scio] ^ (12) 

1 sign(f-£g) I 

IA|-f f-f e/f I A|-f |f — £//[ j 

Focusing on the positive part of the TI spectrum, we see that 
the first term in Jc(££r, 0) contains contributions of the higher 
spin-down band at Sh < t (see Fig. 2). The magnetic field 
shifts the position of this band to higher energies. As a result, 
the tunneling between the leads via virtual states of this band 
is suppressed. The second term includes contribution from 
the low-energy spin-up band. The magnetic field in £h < t 
region lowers the position of the band and increases the tun¬ 
neling probability via its virtual states. This term changes 
sign at Eh = t and becomes negative at large magnetic fields 
Eh > t. The external magnetic field inverts the conduc¬ 
tion band curvature and flips the spin of the conduction band 
within k G [—ko,kQ]. The sign change in the supercurrent 
can be understood by noting the fact that the band-inversion 
permutes the order of tunneling of two electrons through the 
TI and is responsible for reversing the supercurrent direction. 
Further increase of the magnetic field suppresses the tunnel¬ 
ing probability via this band. The negative coupling discussed 
here is similar to the Josephson junction with magnetic impu¬ 
rities studied in Refs. 26-28. 

To summarize, at the charge neutrality point, 6^ = 0, and 
for small eh < t, an increase of the magnetic field increases 
the supercurrent Jc{eh,0)- When eh = t, the critical super¬ 
current changes sign. Further increase of the magnetic field 
decreases the absolute value of the cutTent. The critical cur¬ 
rent as a function of chemical potential and magnetic en¬ 
ergy Eh is shown in Fig. 3. We note that the supercutTent re¬ 
versal can be observed as long as (5/i is smaller than \t — eh\- 
Although we discussed our results within the zero temperature 
limit, a finite temperature only smears the singularities of Eq. 
(11) in the energy interval ^ T around the energy \tztSfi\. 

Having analyzed the supercurrent characteristics, we ap¬ 
ply the condition of weak orbital regime |A| > VHd/2£'^ > 
Eh = v\kh\, where the Fermi velocity in the superconducting 
leads, vp, is usually larger than the Fermi velocity v in the TI. 
Also in order to drive the insulator-semimetal phase transition 


through the external magnetic field, we assume the tunnel¬ 
ing element t is much smaller than |A| so that the condition 
t = Eh i^ satisfied. In this regime, the critical supercurrent 
strongly simplifies and is given by the following expression, 

Jc{.£h, 0 ) = 2et{TT'^VTivsctlf[Q{t - £h) - |^]- ( 13 ) 

We note that the amplitude of Josephson current through the 
semimetallic phase is therefore smaller than that of the insu¬ 
lator phase. It is also worth mentioning that the supercurrent 
reversal discussed above is not affected by the usual Fraun¬ 
hofer response of critical supercurrent to an external magnetic 
field which can be observed in Josephson ‘wide’ weak links 
[29]. 


V. CONCLUSIONS 

Let us now discuss experimental observability and feasi¬ 
bility of the experiment proposed in this paper to reveal the 
supercurrent reversal in 2D TIs. We consider experimen¬ 
tally relevant parameter values for a 2D TI quantum well, 
HgTe/HgCdTe, adopted from the experiment reported in 
Ref. 20. The separation between edge channels was reported 
equal to d = 400nm. Using band structure calculations, the 
hybridization energy associated with this value of d was es¬ 
timated about t = 10/ieV [22]. This hybridization energy is 
of the same order as the proximity induced superconducting 
gap in the edge states which was estimated by Ag < 20/ieV 
[20]. Using the edge mode velocity v = 5 x lO^cm/s we 
find that the Zeeman contribution to Eh can be neglected 
compared to the contribution of Aharonov-Bohm phase gra¬ 
dient which results in a very small critical magnetic field 
id = O.lmT. Therefore, a TI strip of width d = 270nm will 
require an external magnetic field of order H = ImT to show 
the phase transition and supercurrent reversal discussed above 
[22]. These interesting effects can be experimentally verified 
by constructing the configuration shown in Fig. 1 through the 
materials used in Refs. 20 and 21. We are confident that such 
an experiment may pave the way to confirm the edge-mode 
superconductivity. 
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Appendix A: Josephson current 

Here we present more details of the tunneling supercurrent between two superconductors (SC) through a narrow strip of 
2D topological insulator shown in Fig. 1 in the main text. The tunneling between the SC leads and edge states is treated 
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perturbatively. The tunneling Hamiltonian in the interaction representation is given by 


'WTun(T) = 




drdr' [Tr^r'^B a('^> I’)^B,a(T, f') + Tl-.r'„ ('T, l*)®T.a(T, r') + h.C.] 


(Al) 


where q,(t, r) and □.('t, r) are electron spin a =t)i creation operators in the top/bottom (denoted by index j = T/B) 
SC and top/bottom edge (denoted by the same index j = T/B) of the TI, Tj-y is the tunneling matrix element between TI and 
SC which is assumed to be a real quantity, r is the imaginary time, and r = {x, y, z). The imaginary time dependent current 
operator is defined by 


T{t) = ie 


drdr' r)$B,a('r,r') -h.c.]. 


(A2) 


Here the integrals over r and r' run over the edge states and the volume of SC, correspondingly. It is convenient to define a new 
operator as follows. 




Using this new operator, one can rewrite the current operator and the tunneling Hamiltonian, 

I(t) = ielAsir) - AI{t)], 

'hTunir) = ^[Aj(t) + A](r)]. 


(A3) 

(A4) 

(A5) 


To find the Josephson current we take an average of the tunneling current operator in the lowest order of tunneling between 
the superconductor and TI. Therefore, we need to introduce an imaginary time evolution operator 


S{t) = r.rexp| - J 'HTun{r')dT'^. 


(A6) 


At zero bias voltage the single particle current between the SC electrodes vanishes. We also neglect the proximity induced 
superconducting minigap in the TI and set averages q,(t, r)^jy{T', r')) = (Tt-^J ^(r, ^')) = 0 equal to zero. 

The current density in the fourth order of SC-TI tunneling matrix element reads, 

fl/T 


The Green function in the SC leads is given by 


while the edge states of the TI are described by 


dri dT2 dT3 (Tt- Ab (r) Ab (n ) Aj. (r2 ) Aj. (ra )). 

(A7) 

r';r,r') = (7V$j.a('r, r')), 

(A8) 


(A9) 

t'AA') = -{Tr^i,a{TA)^]y{T' ,r')), 

(AlO) 

/;r,r') = ^(r', r')), 

(All) 


Performing a Fourier transformation to the Matsubara frequencies, one finds r, r') = T r, r'), 

where = 7rr(2n + 1), n G Z, and T is the temperature. By substituting the Green functions introduced above into the 
expression (A7), we arrive at the following expression for the Josephson current at zero bias voltage across the junction, 

f ® 

J = — 2eImT E n dv I*!: I’2)^r2,r3G^BT,/37(^n? 1*3? ^4)^r4,r5 

n,a,/3,7,(5 i —1 

^ '^T,75 1 ^5 ; ^6)^r6 ,rj ^TB,i5q( (^n ; ^7 j ^s)^r8 ,ri ■ 


(A12) 
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This expression simplifies in the tunneling barrier model between the SCs and TI. Noting that the tunneling process takes place 
at the SC-TI boundaries we can write for top/bottom edges, = to{x)6{z)d{y =F d/2)S{r — r') and obtain 

f 

J = — 2eImT E n rfajj TIb , a/3 (Wn; a; 1, a;2) f 0 (a^2) Get. /37 (Wn; 2/2 , a/a) to (2:3) 

n,a,^,7,5 i—1 

^ F^^^s{yin]X^,X4)h[x4)GT-R,Sa{'-^n]Xi,Xi)to[xi). (A13) 

Here, we simplify our notation by writing; F^,^g{uJn-,X3,Xi) = (xa, d/2, 0), (a;4, d/2, 0)), FB,a;3(wn; a/i, 2/2) = 

^B,a/3(wn;(a:i,-d/2,0),(a;2,-d/2,0)), and Get,/ 37(w„; (2/2,-d/2,0), ( 2 : 3 , d/2,0)) = GTB,5a(wn; 2/4,a/i) = 

GxE,(5a(tUn: (2/4i d/2, 0), (Xi, d/2, 0)). 


Appendix B: Green fnnction in the topological insulator 


The Green function in the TI is written in the frequency and momentum representations, in the particle-hole, top-bottom edge 
pseudospin, and spin spaces as follows, 

G(a;„;fc) _ 0 V \ _ fGTT{^n',k) GTB{^n]k) 

0 G(a;„;fc)y’ ’ \GBT(w„;fc) GBB{^n',kl) 

where 



G{uJn\k) = 


Gijk) 


f k) k)\ 

k') Gij^^^{ojji] k') J 


(B2) 


These functions satisfy Eq. (3) in the main text. It is apparent that the Green function in the TI is spin-diagonal. In our 
calculations, we need top-bottom edge pseudospin off-diagonal components of the Green function which is defined by 

GTB,n/u(*^»; = GTB,tA(*‘^n; k) = GsT^^/iiiuJn, k) = ^ (B3) 

GTB,n/u(*‘^"; - GTB.tA(*w„; k) = Gbt.tA = (ia;„ - erY - v'^k'^ - F' 
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